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NITRIC OXIDE FORMATION IN GAS TURBINE ENGINES: 

A THEORETICAL AND EXPERIMENTAL STUDY 

by Thomas Mikus , John B. Heywood, and R. Edward Hicks^ 

1. Introduction 

The gas turbine engine is the dominant power plant in current 
aircraft- It has also achieved widespread use in industrial and 
utility applications. It even has the potential to achieve prominence 
as an engine for automobiles (1). In each of these categories, 
however, there is concern about the emission of oxides of nitrogen 
(NO^) from gas turbine engines. For example, in the case of aircraft 
turbines, the predicted growth in air traffic and the trend toward 
larger, higher pressure ratio engines are expected to greatly increase 
N0 x emissions in the vicinity of major airports. The Environmental 
Protection Agency has, as a consequence, proposed NO^ emission 
standards for aircraft gas turbine engines for January 1, 1979 (2). 
These standards will require about a factor of three reduction in 
emissions (in lb^ NO^/lOOO lb^ thrust-hr, averaged over a landing- 
take off cycle) from current jet engine levels. 

It is well known that gas turbine NO^ emissions are most important 

at high engine power settings. Experimentally, it has been found that 

increases in combustor air inlet temperature and pressure, and thus 

engine compressor pressure ratio, increase NO^ emissions. For a 

given engine operating condition, the fuel injection technique, the 

air distribution, the combustor liner pressure drop and the residence 

time all affect NO emission levels, 
x 


Present address. Corporate Research Laboratory, Exxon Research & 
Engineering, Linden, New Jersey 07036. 

^Present address, Chemical Engineering Research Group, Council for 
Scientific and Industrial Research, Pretoria 0001, South Africa. 
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The oxides of nitrogen which leave the engine are almost all 
nitric oxide, NO, with only a few percent nitrogen dioxide, NO 2 . 

It has been shown that the engine exhaust NO concentration depends 
on rate-limited reactions occurring in the burning gases in the high 
temperature zones of the combustor. As additional air is mixed with 
the combustion products to cool the gases before entry to the turbine, 
the NO formation chemistry freezes. After this point, the NO con- 
centration only changes due to dilution. 

In attempts to achieve required reductions in NO^ emissions, several 

new combustor concepts are being developed. A computer program which 

models the NO^ formation process in gas turbine combustors is described 

herein. The model was developed to assist in optimizing the design of 

combustors. It can be used to extrapolate test data from a combustor 

on a test stand to actual engine operating conditions. It can also 

assist in identifying key design parameters which have the greatest 

effect on NO emissions, 
x 

This type of model has been used extensively to predict N0 x emissions 
from conventional combustors. The results of these studies have been 
reviewed recently by Heywood and Mikus (3) , and it was shown that 
relatively simple models which include only the major parameters affect- 
ing NO formation can be used to predict NO^ emissions. 

The design goals of a number of advanced combustors now being 
developed for low NO^ emissions can be summarized as follows. The 
intent is to create a primary combustion zone within the combustor 
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which is fuel lean, where the fuel and air are well mixed together, with 
as short a residence time in the high temperature zone as is consistent 
with the required combustion efficiency. These characteristics, for a 
given engine operating condition, promise minimum N0 x emissions. 

One of the advanced combustor concepts, the NASA swirl-can modular 
combustor (4) , is built upon the use of discrete modules whose character- 
istic size is small compared to overall combustor dimension (5 cm compared 
to 100 cm) . A recirculating wake is produced downstream of each module 
which can be made lean and intensely mixed; its size is similar to the 
module 1 s characteristic size. This is how the design goals mentioned in 
the previous paragraph are approached. 

The scope of this report is to develop a calculational model for 
the formation of NO^ in gas turbine engines. This model is general in 
nature and should be applicable to a range of combustor concepts and 
combustor geometries; however, the examples and discussions are generally 
directed toward aircraft applications. Emphasis is placed on the applica- 
tion of this model to the NASA swirl-can modular combustor, and comparisons 
of modeling prediction with experimental N0 x data are now available. The 
report also makes immediate use of predictions generated by this NC> x 
model. Experimental data are obtained to support assumptions and para- 
meter values used in the model, so that the modeling predictions may 
be used with some confidence. 

The next section describes the particulars of the N0 x model, including 
the "building blocks" which can be assembled to generate NC> x predictions 
for a wide variety of combustors. The third section assembles these 
blocks and compares the resulting predictions with available N0 x data 
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for different combustors, noting the influence of the modeling parameters. 
Section 4 is an analysis of Section 3 ? s predictions relevant to the 
NASA swirl-can combustor, demonstrating that much can be learned about 
a combustor’s operation simply by matching its NO^ emissions data. The 
fifth section describes parametric modeling studies which provide examples 
of the further usefulness of the NO^ model in extrapolating test data 
and in designing for low NO^; further insight is also obtained into the 
influence of the modeling parameters. 

Sections 6 through 10 are directed toward the verification of the 
assumptions and parameter values which went into the NO^ modeling of the 
NASA swirl-can combustor. The sixth section presents a calculation of 
the streamlines downstream of a module in order to get an overall view 
of the flow in this combustor. Section 7 describes calculations of the 
fuel droplet formation and evaporation processes which occur within that 
flow field. Flow visualization experiments which document the turbulent 
recirculation and fuel droplet processes downstream of the module are 
exhibited in Section 8. The ninth section is a presentation of non- 
burning tracer experiments which quantitatively follow the streamlines 
from the swirler jets downstream, as they entrain adjacent gases. 
Experiments conducted under burning conditions are presented in Section 10 
they map out the flow field in terms of both local equivalence ratio and 
local nonuniformity. A summary and conclusions are set forth in the final 


section 
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2 . Model Description 

2.1 Chemical kinetics of NO formation 

x 

At the temperatures , pressures, and equivalence ratios found 
in gas turbine combustors, NO formation may be calculated from the 
set of six chemical reactions used by Lavoie et al (5) . The three 
most important reactions are a modification of the Zeldovich mechanism: 


N 2 + 0 t N + NO 

(i) 

NO + 0 t- N + 0 2 

(2) 



N + OH NO + H 

(3) 


Reactions involving N 2 O as an intermediary may become important at 
temperatures below 2100 K, especially if concentrations are substanti- 
ally greater than equilibrium values. Three such reactions which 
were included in the calculational scheme are: 


H + NjOt 

n 2 

+ OH 

(4) 

N 2 + 0 2 t 

0 

+ n 2 o 

(5) 

NO + NO t 

0 

+ n 2 o 

(6) 


A reaction cited by LeTrung et al (6) which came under consideration 
too late to be used in many of the calculations is: 


N 2 0 + M£N 2 +0 + M 


(7) 
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For typical conditions of interest the above reaction is more 
important than reactions (4) and (5), but as will be shown below, 
all of the N 2 0 reactions may generally be neglected with little ef- 
fect on the NO formation rate. Rate constants taken from Baulch 
x 

et al. (7) are presented in Table 1. 

The reaction scheme may be written as rate equations for NO, N, 
and N 2 0 (if included), but first the method for calculating the con- 
centrations of 0, 0 2 , OH, and H must be resolved. The N 2 concentration 
is essentially constant. The appropriate assumptions depend on the 
pressure, temperature and atomic composition of the system being 
analyzed. These assumptions should be suitable for typical gas turbine 
combustor conditions, i.e., pressures - 5-30 atm, peak temperatures 
- 2500 K, residence times at these peak conditions = 1-3 msec, 

equivalence ratio for hydrocarbon-air combustion between about 0.6 
and 1.3. 

A model for the flame front in the combustor must be formalized 
before appropriate assumptions can be made. Heywood (8) and Westen- 
berg (9) have justified the assumption that under appropriate condi- 
tions, the 0, 0 2 , OH and H concentrations are the equilibrium values. 
The argument can be summarized as follows. 

In a premixed one-dimensional flame, there is a thin initial 
reaction zone where radicals are generated and react with the hydro- 
carbon fuel. Most of the temperature rise occurs in this zone. This 
region is followed by a thicker zone where CO is oxidized to CO 2 , and 



7 


TABLE 1 

Sate Constants for the NO^ Kinetic Reaction Scheme 


Reaction , i 
1 
2 

3 

4 

5 

6 

7* 


Rate const., k^, cm 3 gmole 3 sec 
7.6 x 10 13 exp (-38,000/T) 

1.5 x 10®T exp (-19,500/T) 

4.1 x 10 13 

7.6 x 10 13 exp (- 7,600/T) 

6.3 x 10 13 exp (-55,200/T) 

1.3 x 10 12 exp (-32,100/T) 

5.0 x 10 14 exp (-29 ,000/T) 

where T is temperature in K 


for M * A ; no recommendation for M * N 2 
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radicals present in greater than equilibrium concentrations recombine** 
through trimolecular reactions. If the NO formed in this non- 
equilibrium region is much less than the amount formed downstream 
of the flame, then the 0, 0 2 , OH, and H equilibrium assumption is 
a reasonable approximation. There have recently been studies of 
CO equilibration in the flame zone at pressures, temperatures and 
equivalence ratios typical of a combustor primary zone (10). These 
studies show that equilibration times are less than 2 x 10 4 sec 
(for equivalence ratios between 0.6 and 1.3) which is much less than 
typical primary zone residence times, and that the effect of radical 
concentrations above equilibrium values within the hydrocarbon reaction 
zone will not significantly increase NO emmissions. The amount of 
NO formed in the reaction zone of the flame is small compared with the 
amount formed downstream of the flame at gas turbine combustor 
conditions. 

Following Lavoie et al (5) and assuming equilibrium concentrations 
for N 2 ,0, 0 2 ,0H, and H, the rate equations for NO, N, and N 2 0 can be 
written as 


4— = -a (BRi + R 2 + Ra + 


p = -a (BRi + R 2 + R 3 + 2aR 6 ) + Rj +B (R 2 + R 3 ) + 2yR 6 (8) 

at P 


_ d [N] 
P dt p 


= — B (OiRi + R 2 + R 3 ) + Ri + 0t (R 2 + R 3 ) 


(9) 


d t M 2°1 

dt p 


= — Y (R 4 + R 5 + R6+ R 7 ) + R 4 + R 5 + a R 6 + R 7 


( 10 ) 
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where p is the gas density; a = [NO] /[NO] , g = [N]/[N] * and y = [N 2 0]/ 

e 6 

[N20]^ are, respectively, the concentrations of NO, N, and N 2 0 divided 
by their equilibrium concentrations; and R^ is the "one-way lf equili- 
brium rate of the ith reaction, e.g., Ri = ki[N 2 ] [0] . Rough estimates 

e e 

of the magnitude of the terms in these equations indicate that the 
relaxation times for the [N] and [N 2 0] equations are several orders 
of magnitude shorter than that for the [NO] equation. It is therefore 
an excellent approximation to assume steady-state concentrations for 
N and N 2 0 and set d([N]/p)/dt and d([N 2 0]/p)dt = 0. Equations (9) 
and (10) can then be used to eliminate g and y from equation (8) ; 
rearranging in terms of the mass fraction of nitric oxide, {NO}, the 
rate equation becomes 


d{NQ} 

dt 


%0 

p 


(1-a 2 ) ( 


Ri 


1+aKi 


Re 


i+k 2 


( 11 ) 


Where {NO} is the mass fraction of NO, is the molecular weight of 

NO, Ki = Ri/(R 2 + R 3 ) and K 2 = Re/(R 4 + R 5 + R 7 ) - The first term in 
the large bracket on the right-hand side of the above equation is the 
result of reactions (l)-(3), the extended Zeldovich chain mechanism. 

The second term is the result of reactions (4)-(7) and involves N 2 0 as 
an intermediary. In general, Ki , K 2 , and a are of order unity or less, 
and the relative importance of the two mechanisms just described is 
determined by the ratio of Ri to R6 • This is much larger than unity 
except at low temperatures (less than 2100 K) . 
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Typical values for the quantities found in equation (11) are given 

in Table 2. Many of the calculations presented herein have been made 

with R^O; however, reactions (4)- (7) usually contribute relatively 

little to N0^ formation, as can be seen from an examination of equation 

(11) in light of the values presented in Table 2. For a local mixture 

richer than (J> ^ 0.7, the R,, term is much smaller than the R. term, barring 

o 1 

rare situations where NO is far above equilibrium levels (ot»l). For 
a local mixture leaner than <J> ^ 0.6, the R^ term becomes comparable to 
the R^ term, but both are insignificant relative to richer mixtures. In 
addition, for mixtures which are too lean, the kinetic scheme is invalid 
due to the combustion reaction time constants becoming comparable to the 
residence time. Thus, except for cases where essentially all (on a mole- 
cular scale) of the burned gases are quite lean, equation (11) may be 
simplified to 



Equations (11) and (12) show that the NO formation rate depends only 
on the pressure, temperature, atomic composition and NO concentration 
downstream of the thin flame zone. NO profiles calculated with equation 
(11) as a function of time downstream of a one-dimensional flame for 
a stoichiometric mixture are shown in Fig. 1. Since primary zone resi- 
dence times are of order or less than 3 msec and peak temperatures ^ 

2500 K, one would not expect NO concentrations to reach equilibrium 
levels within the high temperature zone of a gas turbine combustor. 

This argument can be made more precise. For a << 1, equation (11) 
can be rearranged as 
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TABLE 2 


Typical Values for Quantities Entering 

into Rate Equation (11) 

for NO 

inlet air temperature 

590 K 



pressure 

5.8 atm 



fuel 

kerosene 




30 g gmole 

■1 


equivalence ratio 

0.8 

0.5 


burned gas temperature, K 

2250 

1730 


-3 

P» g cm 

9.1 x 10" 4 

1.2 x 

10~ 3 

_ , -3 -1 

R. , gmole cm sec 

7.5 x 10 -7 

2.3 x 

10" 10 

K 1 

0.41 

0.20 


n . -3 -1 

R^, gmole cm sec 

3.0 x 10 -8 

1.1 X 

io - 10 

K 2 

0.029 

0.008 


K 2 for R ? = 0 

0.21 

0.63 


dim 8 -1 

dt ’ S£C 

5.1 x 10" 2 

1.7 x 

IQ " 5 
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1 d{NO} = 2M N0 R 6 . _ 1 

{NO} dt pTno} ^1 1 + K„ ; T„ 

e e / no 


(13) 


{NO}^ 1 d{NO}/dt has the units sec ^ and can be thought of as the recipro- 
cal of the characteristic time, t N q» °f the NO formation chemistry. 
Generally, this equation may be simplified to 


1 d{NO}% 2M N0 R 1 a, 1 

{NO} dt plNO} T„_ 
e e NO 


(14) 


as was argued for equation (12) . 

Since the gas turbine combustor is essentially adiabatic, the burned 
gas temperature and concentrations of major species are determined by 
the combustor inlet air conditions, local equivalence ratio, and fuel 
properties. In a simple-cycle gas turbine engine, combustor inlet 
conditions are determined by engine inlet conditions, compressor pres- 
sure ratio and efficiency, and pressure losses. Values of t. 7 _ from 

NO 

equation (13) are plotted in Fig. 2 as a function of equivalence ratio 
for different engine compressor pressure ratios, assuming 1 atm., 288 K 
inlet air, 90% polytropic compressor efficiency, and 3% total pressure 
loss across the combustor liner (see Table 3). 

The kinetic scheme for NO formation presented in this subsection 
will be tested in Subsection 3.1. 


2.2 Effects of fuel-air nonuniformities; a primary zone model 

To apply the above NO kinetic scheme to a heterogeneous combustion 

zone, this description of a premixed one-dimensional flame must now be 
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TABLE 3 

Inlet Temperatures and Pressures Appropriate- to the NO^ Time Constants 
and Formation Rates of Figures 2 through 6 


Compressor 

Pressure 

Ratio 

Combustor 
Inlet Air 
Temp, K (°F) 

Combustor 
Pressure 
atm (psia) 

8:1 

555 (538) 

7.6 (112) 

16:1 

685 (774) 

15.2 (223) 

24:1 

772 (933) 

22.8 (335) 

32:1 

840 (1052) 

30.4 (447) 



MEAN PRIMARY ZONE EQUIVALENCE RATIO, </> p 

Rate of Increase in NO Emission Index versus 
Mean Equivalence Ratio for Various Values of 
Nonuniformity Parameter, 8:1 Engine Pressure 
Ratio; see Table 3 for inlet conditions. 


Figure 3. 
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related to the flame structure in a gas turbine combustor. Fuel as a 
liquid spray and air enter the primary zone separately in approximately 
stoichiometric proportions. The fuel droplets vaporize as they move 
relative to the primary zone gases leaving wakes rich in fuel vapor, which 
then mix with air and already burned gas, and bum. Mixing of fuel and 
air to a uniform equivalence ratio would not be expected. As the fuel 
droplets vaporize, fuel vapor and air mix and form turbulent "eddies" or 
pockets of combustible mixture with a wide range of fuel-air ratios. 

High speed movies of liquid-fueled combustors and burners suggest 
that the thin flame front model described in Subsection 2.1 is still 
applicable. In such combustors low luminosity blue flame fronts propa- 
gate through some of the fuel/air eddies, while soot-forming, luminous, 
yellow flame fronts envelope and consume other fuel/air eddies. Such 
an unsteady turbulent flame structure would be expected in the vigorously 
stirred primary zone where fuel and air enter separately. As a conse- 
quence, burned gas eddies are produced with a distribution in fuel-air 
ratio about the mean primary zone value. This distribution will change 
with time as these eddies mix with each other, and with dilution air, 
through the action of turbulence and molecular diffusion. 

Details of the fuel introduction method and air flow pattern de- 
termine the initial droplet size and spatial distribution. The air 
flow pattern further dictates the process of droplet motion and evapora- 
tion as well as the subsequent mixing intensity. Thus, fuel introduction 
and air flow affect both the initial spread of the burned gas fuel-air ratio 
distribution and its change with time. This physical description of the 
combustion process is the basis for the model which will now be described. 
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The importance of nonuniformities can be illustrated as follows. 
Since fuel and air are not uniformly mixed, different burned gas eddies 
Will have different average fuel fractions during their residence time 
in the combustor primary zone. The use of fuel fraction F (fuel mass 
divided by total mass) in a mass-based distribution function simplifies 
the mathematics, compared to using fuel-air ratio or equivalence ratio. 
Assuming a Gaussian distribution function for fuel fraction about the 
mean fuel fraction F, 

f (F) = (l/2TTCf 2 ) 1/2 exp{-(F - F) 2 /2a 2 } (15) 

then the fraction (by mass) of the burned gas flow with fuel fraction be- 
tween F and F + dF is f (F) dF. a is the standard deviation of the dis- 
tribution and a mixing parameter 

s = O/F (16) 


can be used as a measure of nonuniformity. 

Figures 3 through 6 show the effect of nonuniformities on NO forma- 
tion rates (plotted as dE^ /dt) for four different engine compressor 

pressure ratios; E is a global emission index expressing oxides 

NO. 2 

of nitrogen as equivalent per 1000 kg of fuel burned. The 

simplifying assumption has been made that NO is far below equilibrium 
(a « 1). Equation (11) has been used to calculate the perfectly 

mixed case (s * 0) . For nonuniform mixtures (s > 0) values of dE^^/dt 
obtained from equation (11) for a range of cj> are weighted by equation 


I 
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(15) to determine an average NO formation rate at each average equiva- 
lence ratio. For stoichiometric mixtures, the mean NO formation rate 
is substantially reduced as the primary zone becomes less uniform. For 
rich 1.2) and lean (<t>^< 0.7) the mean NO formation rate is substan- 

tially increased. It has been shown that values of s between 0.3 and 0.7 

are required to match predicted N0 x emissions to measurements in conven- 
(3) 

tional combustors . 

Approximate estimates of combustor N0 x emissions can be made using 
Figures 3 through 6. For example. Figure 5 is appropriate for combus- 
tors used In state-of-the art engines in the 40,000 + lb^ thrust class 
at take-off conditions. Assuming an s = 0.4 is typical for the combus- 
tor, dE„_ /dt is almost constant between 0.6 < TfT < 1.0. Assuming a resi- 

„ NO 2 p 

dence time for the burned gas within this equivalence ratio range of 1 
to 2 msec, and an average value for dE^ /dt of 24 kg NO^/IOOO kg 
fuel msec, gives an E^ between 24 and 48 kg NO z/1000 kg fuel. This 
spans the measured range of 30-40 kg N0 2 /1000 kg fuel. 

These graphs also indicate desirable combustor operating conditions 
for low N0 x emissions, dE^ /dt decreases below levels found in 
conventional combustors only for (f>^ ^ 0.7, and then only if the uniformity 
of the fuel-air mixture is significantly improved. Thus a lean, well 
mixed, short residence time primary combustion zone is required, with 
the most significant gains coming from better fuel-air mixing. 

For detailed N0 x calculations a primary zone treatment has been 
implemented using the approaches set forth in this subsection. The 
primary zone is modeled as a partially stirred (or macro-mixed) reactor. 

It is assumed that the fuel and air which enter the primary zone are 
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Figure 4. Rate of Increase in N0 X Emission Index 
versus Mean Equivalence Ratio for 
Various Values of Nonuniformity Parameter 
16:1 Engine Pressure Ratio; see Table 3 
for inlet conditions. 
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Figure 5. Rate of Increase in N0 X Emission Index 
versus Mean Equivalence Ratio for 
Various Values of Nonuniformity Parameter. 
24:1 Engine Pressure Ratio; see Table 3 
for inlet conditions. 
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Figure 6. Rate of Increase in N0 X Emission Index 
versus Mean Equivalence Ratio for 
Various Values of Nonuniformity Para- 
meter, 32:1 Engine Pressure Ratio; 
note N0 X scale shift from previous 
figures; see Table 3 for inlet 
conditions . 

TABLE 4 

Typical Data Appropriate to the 
Premixed N0 X Calculations of Figure 7 


unburned mixture temperature 

pressure 

fuel 

equivalence ratio 0.7 

burned gas temperature 2380 K 


800 K 
5.5 atm 


propane 

1.0 

2670 K 


residence time 


3.6 msec 


3.1 msec 
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rapidly dispersed throughout the zone to form discrete fluid eddies of 
a scale small compared with combustor dimensions. 

For a stirred reactor, the fraction (by mass) of the flow with 
residence time between t and t + dt is ip(t)dt where: 


t) = (1/ T ) exp(-t/x p ) 


(17) 


Tp is the mean residence time for the primary zone, which can be calcu- 
lated as 


T 


P 



(18) 


where is the primary zone volume, m^ in the mass flow rate through 
the zone, and v is the mean specific volume within the zone. Since 
all of the fuel must pass through the primary zone. 


m 

P 


F 

o o 


F 


P 


(19) 


where m is the overall combustor mass flow and F /F is the ratio 
° p o 

of primary zone to overall fuel fraction. By definition 

v * | v(F) f(F) dF (20) 

0 

where f(F) comes from substituting F = F^ and s » s^, the primary zone 
nonuniformity parameter, into equations (16) and (15). Thus, three 
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parameters are required to model the primary zone. V^, F^/F^ and 
were chosen because they are expected to be weak functions of combustor 
operating conditions. 

NO^ predictions for the primary zone are then calculated as 


{NO} *= f f {NO}(F,t) dt f (F) dF 

P J Q * 0 


( 21 ) 


where {NO}(F,t) is obtained by integrating equation (11) starting 
from {NO} = 0 at t = 0. {NO} might also be started from a level which 
approximated "flame-front NO", but this level is negligible for the 
cases presented herein. The integral on t has been derived semi-analy- 
tically. The integral on F is performed by creating a statistical en- 
semble of equal mass "fluid elements" whose fuel fractions represent 
f(F), and then averaging over these elements. 


2.3 A multi-zone flow model with Monte Carlo mixing 

In general, NC> x production is not confined to the primary zone. 
Further downstream where additional air is added, it is necesary to 
model one or more zones where mixing processes are proceeding con- 
currently with N0 x formation. A statistical method of following a 

chemical reaction (such as nitric oxide formation) through the course 

( 12 ) 

of a mixing process was described by Corrsin , and has been adapted 

(13) 

more recently by Flagan and Appleton 

With this method the flow through a zone is represented by a sta- 
tistical ensemble of fluid elements of equal mass. At any point in 
time each element has its own local fuel fraction. The number of 
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elements should be large enough to fairly characterize the distribution 
of fuel fractions. This number was between 250 and 1000 for the calcu- 
lations presented herein. The ensemble of fluid elements experiences 
mixing interactions which represent the dilution and mixing process. 

These interactions occur at time intervals of 

It = 1/BN (22) 

where B is the current value of the mixing intensity and N is the total 
number of fluid elements in the statistical ensemble. Each mixing in- 
teraction is calculated as follows: (i) Two different fluid elements are 
chosen at random. (ii) For each fluid element the nitric oxide formation 
parameters corresponding to its local fuel fraction are obtained by inter- 
polation from a table of parameters appropriate to the combustor inlet 
air pressure and temperature and the fuel properties. These parameters are 
the specific volume of the burnt gases, v, and a set of parameters 
required for the use of equation (11). (iii) Starting with the nitric 
oxide mass fraction each element contained immediately after its previous 
interaction, equation (11) is integrated semi-analytically up to the 
current time. If this is an element's first interaction, the integration 
starts with the nitric oxide mass fraction contained when it entered the 
zone, (iv) These integrated nitric oxide mass fractions are averaged to 
represent complete mixing of the two interacting elements. (v) Both 
elements are assigned a new fuel fraction equal to the average of the 
two prior to mixing. 

For each time increment there is a corresponding volume increment 
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which represents the progress of the flow through the zone's volume 

AV = m At v (23) 

where m is the zone's mass flow rate and v is the current average speci- 
fic volume. When the entire zone volume has been swept out by these 
incremental volume flows, the local nitric oxide mass fractions for 
all of the fluid elements are updated and averaged. This average re- 
presents the mean nitric oxide mass fraction at the zone's exit. 

Within the framework of this model is the freedom to specify the 
mixing intensity 3 as a function of time. This will be utilized in 
the modeling calculations presented in Subsection 3,2. 

The dilution zone model is easily linked into a computation. If 
the preceding zone is a primary zone model as described in Subsection 
2.2, there will be an ensemble of fluid elements corresponding to the 
flow which passes through the primary zone. Each element has its own 
local fuel fraction and current nitric oxide mass fraction. To form 
the ensemble of fluid elements for the dilution zone is simply a matter 
of taking the primary zone ensemble and adding a number of new elements 
which corresponds to the mass flow rate of additional air which is 
entering the dilution zone. These new elements will of course have 
fuel fractions and nitric oxide mass fractions equal to zero, since 
they represent pure air. 
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3- Comparisons of NO^ Predictions with Available Experimental Data 

3.1 A "premixed" tubular burner 

The kinetic scheme for nitric oxide formation described in 
Subsection 2.1 is best tested by comparing predicted NO levels with 
measurements from a premixed burner. The fuel and air should be 
completely premixed (s=0) , and the flow should be adiabatic, with a 
well-defined residence time - a one-dimensional flow system, for 
example. Anderson (14,15) presented NO^ measurements from a burner 
which approached these requirements. It was a tubular burner 4 
inches in diameter, consisting of a propane inj ector/mixer , 20 to 
25 diameters of mixing region, a plate with cylindrical passages for 
flameholding and up to 4.5 diameters of reaction zone; a cooled gas 
sample probe was located on the axis of the burner. 

For this geometry a plug flow assumption should be quite accurate. 
The radiative heat transfer from the burner centerline was estimated, 
and found to decrease the burned gas temperature by less than 15 K; 
therefore, the gas flowing along the centerline was assumed to be 
adiabatic. The triangles in Figure 7 are the NO^ emission index data 
at various equivalence ratios presented by Anderson (15), measured 
30 cm downstream of the flameholder, for the case of 5.5 atm, 800 K 
inlet mixture temperature, and 25 m/sec reference velocity. The curve 
labeled A shows the NO^ model predictions for the completely uniform 
case; some data relevant to these predictions are given in Table 4. 
Curve A is generally in good agreement with the data over a wide 
range of emission index, especially in the main region of interest 
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10 < < 40, The NO^ predictions are not expected to be accurate 

at equivalence ratios lower than about 0.5 f where temperatures are 
lower than 2100 K (see Subsection 2,1). At equivalence ratios near 
unity Anderson states that his data may be in error due to catalytic 
reduction of NO in his stainless steel probe; this would explain some 
of the discrepancy between predictions and data at the higher 
equivalence ratios. 

As seen in Subsection 2,2, a less peaked prediction curve may 
be generated by assuming a certain degree of nonuniformity (s>0) . 
Curve B of Figure 7 shows NO^ predictions for s = 0.1. This curve 
follows the general trend of Anderson’s data significantly better 
than curve A. It should be noted that these data were obtained 
downstream of a water-cooled flameholder, which would lower burned 
gas temperatures; the predictions assume adiabatic conditions. If 
curve B were corrected for heat transfer to the flameholder, it 
would agree with the data even more closely. 

Figure 8 presents comparisons for lower inlet mixture 
temperatures. The triangles in this figure are again data from 
Anderson (15); conditions were the same as for the data of Figure 7, 
except that the inlet mixture temperature was 600 K. The solid 
curve labeled A again shows NO^ model predictions for the completely 
uniform case, while solid curve B is for s « 0.15. The agreement 
between the predictions and the data is not as good as for the higher 
inlet mixture temperature case, but as in Figure 7 agreement is 
better for the nonuniform predictions. Burned gas temperatures will 
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be below 2100 K for equivalence ratios less than about 0.6 in this case, 
so that NO^ predictions are not expected to be accurate in that 
range. As in Figure 7, catalysis in the sample probe may have 
depressed the NO^ data at the richer equivalence ratios, and a 
flameholder heat transfer correction would lower the predictions; 
these factors would improve the agreement. 

The circles in Figure 8 are data from an earlier Anderson paper 
(14). These data were measured 46 cm downstream of the flameholder, 
for the case of 5.5 atm, 590 K inlet mixture temperature, and 
23 m/sec reference velocity, roughly comparable to the other data 
in Figure 8. However, the data circles were obtained with a 
different propane injector /mixer design, and an uncooled flameholder 
only one-third as thick. The dashed curves labeled C and D are the 
NO^ emission predictions for s=0 and s=0.1, respectively. These 
predictions follow the data quite nicely up to an equivalence ratio of 
about 0.85; above this level there is some discrepancy. Sample probe 
catalytic effects could explain some of this. Generally, this set 
of data seems more uniform, i.e. is better matched by NC> x predictions 
with a low s value, probably due to differences in the propane 
injector/mixer. Also, heat transfer to the flameholder would be less 
significant than for the other data, so that the NO^ predictions need 
not be adjusted downward. 

There are two conclusions which may be drawn from a comparison 


of the predictions and data presented in Figures 7 and 8. (i) The 

nitric oxide kinetic scheme agrees reasonably well with measurements 
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from a simple burner. (ii) Although it was much more uniform than 
state-of-the-art gas turbine combustors, the "premixed" burner 
described above probably did not produce a completely uniform fuel-air 
mixture. 

For a stoichiometric mixture, a calculation of the average 
mixing intensity required to produce a nonuniformity of s=0.15 
at the flameholder (starting from pure fuel and air at the fuel 
injector) yields the value 3=70 sec Compare this number with a 
value of approximately 20 sec ^ for pipe-flow turbulence under these 
conditions. The additional mixing could easily be accounted for by 
considering the initial mixing produced by the propane inj ec tor /mixer . 
Thus, the fuel-air mixture in this burner should not be expected to 
be completely uniform. 

3,2 The NASA swirl-can modular combustor 

A more complete evaluation of the NO^ model described in Section 2 
requires the examination of a combustor with a flow field more 
complicated than the burner described in the previous subsection. 

The NASA modular swirl can combustor is an unconventional annular 
design (4). Fuel is introduced into 120 swirl cans arranged in three 
concentric rings at the same axial station. Most of the air flows 
through these modules and in between them at the plane of the 
"blockage plates" which are the maximum axial cross-sections of the 
swirl cans. Only about 6% of the total air flow is used as liner 
cooling air and does not pass through or between the swirl cans. 


i 
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Figure 9 shows the construction of a swirl can. Air flows from the 
left, entering the can and mixing with fuel, which is sprayed onto 
the swirler plate from a fuel tube inside the can. The fuel-air 
mixture passes through the swirler and burns in the wake of the 
module. Some of the air which has passed around the module is 
entrained into the recirculating wake flow. Cold-flow tests confirm 
that this recirculating wake extends approximately one blockage-plate 
diameter downstream. The air which flows around the module but is 
not entrained into this "primary zone" will mix with the combustion 
products from the primary zone further downstream. Figure 10 is a 
tangential view of the combustor showing the swirl can modules 
mounted at the entrance plane of the annular liner. 

At low fuel-air ratios the combustion occurs within small 
primary zones in the wake of each swirl can. As the overall fuel-air 
ratio approaches stoichiometric (0.069), the primary zones coalesce 
into a single "ring of fire." These two extremes of fuel-air ratio 
serve as limiting cases which are the bases for the first two sets 
of NO^ predictions described in this subsection. These limiting 
cases allow the NC> x calculation to be made with a single "primary" 
zone as described in Subsection 2.2, Later in this subsection the 
dilution zone presented in Subsection 2.3 will be added to make a 
more complete model of NO^ formation in this NASA combustor. 

In the limit of low fuel-air ratios, so little fuel is being 
added to the flow through the swirl cans that combustion in the 
module wakes occurs at fuel-lean conditions. Any subsequent mixing 
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with the air which passes between the modules but is not entrained 
into the Immediate wakes will make the combustion products even leaner, 
rapidly decreasing the rate of nitric oxide formation. This effect 
can be seen in Figures 3 through 6. Regardless of which curve is 
applicable to the case of interest, at sufficiently low equivalence 
ratios the nitric oxide formation rates fall off sharply as the 
equivalence ratio decreases. Therefore, the complete nitric oxide 
formation process can be modeled in this limit as the products of a 
stirred-reactor module-wake primary zone being suddenly quenched and 
then simply diluted to the overall fuel-air ratio. This model 
corresponds to the separate module wake combustion which is observed 
at low fuel-air ratios. 

An uncertainty encountered in analyzing the NASA combustor is that 

contrary to more conventional designs, the volume of the primary zone 

and the mass flow of air entering this zone are not well defined by 

the liner geometry. The fluid mechanics of entrainment into a 

reacting, swirling, recirculating flow are not tractable enough to be 

used in the relatively simple nitric oxide formation model being 

proposed here. The following assumptions are made to obtain simple 

models: (i) The volume of the primary zone in the wake of each 

3 

swirl can module is taken to be 8.7 in , which is equal to the area 
of an averaged blockage plate multiplied by the plate ? s diameter. This 
aspect ratio of approximately one is verified by cold-flow tests. If 
the flow patterns in the combustor do not change much as the volume 
flow rate and combustion intensity are varied, the volume of the 
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primary zones will remain nearly constant. Of course, as the fuel- 
air ratio approaches stoichiometric and the primary zones coalesce, 
this assumption becomes invalid, (ii) The average fuel fraction in 
the primary zones is proportional to the combustor’s overall fuel 
fraction. Once again, if the combustor’s flow patterns remain 
relatively fixed, this assumption will be approximately true, 

A few comments on the above modeling assumptions are in order, 
(i) Primary zone residence times are proportional to primary zone 
volume. Nitric oxide production is found to be almost linear with 
residence time for the region of interest, so that changing the 
primary zone volume will approximate sliding the nitric oxide 
prediction curves vertically on a log scale. (ii) Changing the 
proportionality factor between the primary and overall fuel fractions 
is equivalent to changing the fraction of the overall air flow which 
enters the primary zone. This changes both the primary zone 
residence time and the dilution factor on exit from the primary zone. 
The approximate linearity of nitric oxide with residence time 
mentioned above results in a cancellation of the residence time and 
dilution factor effects. Thus, changing the proportionality factor 
between the primary and overall fuel fractions merely slides the 
nitric oxide predictions horizontally on a log scale of overall 
fuel-air ratio. 

In the limit of overall fuel-air ratios near stoichiometric, so 
much fuel is being added to the flow through the swirl cans that 
combustion conditions in the module wakes are much richer than 
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stoichiometric. As mixing occurs with the air which passes between the 
modules but is not entrained into the immediate wakes, the rate of 
nitric oxide formation will increase markedly. This effect is shown 
in Figures 3 through 6. Regardless of which curve is applicable to 
the case of interest, at equivalence ratios greater than unity, the 
nitric oxide formation rates are significantly lower than the rates 
prevailing after dilution to an equivalence ratio near unity. It is 
also found that residence times calculated for the module wakes are 
several times shorter than the overall residence times calculated for 
the combustor. For the above reasons, the nitric oxide formed in 
the module wakes can be neglected compared to that formed after 
dilution to the overall fuel-air ratio. Therefore, the complete 
nitric oxide formation process can be modeled in this limit as the 
whole combustor volume operating as a stirred reactor at the overall 
fuel-air ratio. This model corresponds to the "ring of fire" type of 
combustion observed at fuel-air ratios near stoichiometric. 

Since the entire combustor liner annulus is serving as the 
combustion zone of interest, the volume of the zone as well as the 

mass flow rates through it are well defined by the liner geometry. 

3 

The volume of the annulus is approximately 4.1 ft . The average fuel 
fraction is about 1.03 times the overall to allow for one half of 
the liner cooling air to participate in the reactor on the average. 

Predictions from the two limiting case models have been 
published previously by Heywood and Mikus (3) • They serve as an 
important backdrop for the more complete model which is implemented 



34 


later in this subsection. Predictions are expressed as nitric oxide 
mass fraction {NO}. These may be converted to the emission index 
by the formula 


"NO, 


= 1.53 


{NO}(l+f) 

f 


x 10“ 




1.53 


{NO} 


10 “ 


(24) 


where f is the overall fuel-air ratio. Calculations have been done 
for both the residence time distribution of equation (17) and for a 
uniform residence time; these predictions differ by only a few percent. 
The latter method has been used for the predictions presented in 
this subsection. 

Figures 11 through 13, and 14, show these limiting-case predictions 
as dashed lines, for two different inlet air temperatures. Data 
points shown in these figures were derived from NASA data (4). The 
data points shown as circles are measurements made at a combustor 
inlet pressure of 5 to 6 atm, with a combustor air flow of 85 to 110 
lbm/sec. The data points shown as diamonds are known to be measured 
at 6 atm and 110 lbm/sec, which were the values used in the prediction 
computations. The solid curves in these figures include the dilution 
zone model, and will be described later in this subsection. 

The dashed curves drawn for fuel-air ratios between 0.01 and 0.03 
are calculated with the low fuel-air ratio model, using a primary- zone 
fuel fraction 2.5 times the overall value. As mentioned earlier, a 
change in this parameter will simply proportion out the prediction 
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FUEL - AIR RATIO 


Figure 11. Comparison of N0 X Predictions with 
NASA Swirl-Can Data (4) for Various 
Values of a Constant Mixing 
Intensity (solid curves); fixed 
primary zone parameters. 
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Figure 12. Comparison of N0 X Predictions with 
NASA Swirl-Can Data (4) for Various 
Values of Mixing Intensity Decay 
Time (solid curves); fixed primary 
zone and initial intensity. 
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curves horizontally. Also note that a change in primary zone volume 
simply shifts the curves vertically. The mixing parameter values 
0.33 and 0.42 which are illustrated produce approximately the proper 
slope to match the data, especially for the higher combustor inlet 
temperature case. Figures 11 through 13. 

The dashed curves drawn for fuel-air ratios greater than 0.03 
are calculated using the model appropriate to fuel-air ratios 
comparable with stoichiometric. Here the primary zone volume and 
fuel— air ratios are well defined, so that the prediction curves are 
not shiftable without making basic changes in the model. 

Having matched the emissions data with the limiting-case models, 
the following model parameters can be assigned approximate values for 

o 

this combustor: (i) primary zone volume V £ 8.7 in. (ii) primary 

P 

to overall fuel fraction factor /F q ^2.5 and (iii) primary zone 
mixing parameter s^ % 0.4. It should be noted that these parameters 
were only evaluated approximately, and that they by no means represent 
the best possible match with the experimental data. However, these 
parameter values are adequate for exploring the potential of the 
dilution zone model to predict nitric oxide formation over a wider 
range of combustor fuel-air ratios. 

The dilution zone calculation is in fact an extension of the lean 
end primary zone calculations. As described at the end of Subsection 
2.3, the dilution zone is linked in by adding additional fluid 
elements for the dilution air. In this case these new elements 
represent the air flow which has passed between the modules, but 
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was not entrained into the recirculating wake primary zone. For a 

primary zone enrichment factor F /F =2.5, only about 40% of the 

p o 

air flow passes through the primary zone. So for a dilution zone 
ensemble of 750 elements, about 300 are carried over from the 
primary zone ensemble, and 450 will start out as additional air 
elements. Starting from this condition, the mixing and NO 
calculations are carried forward as described in Subsection 2.3. 

Figure 11 presents the mixing model predictions for a range of 
values of the mixing parameter 3, assuming that in each case 3 is 
constant in time. At a given fuel-air ratio, the behavior of the 
nitric oxide predictions as a function of 3 can be understood by 
examining the history of the nitric oxide formation rate as the 
flow from the primary zone is mixed with the dilution air. For 
example: (i) At overall fuel-air ratios less than about 0.02, the 

primary zone is leaner than the value for peak nitric oxide formation. 
As the flow from the primary zone is diluted, the nitric oxide 
formation rate decreases sharply, becoming negligible long before the 
overall fuel-air ratio is reached. Therefore, a high value of 3 
representing rapid dilution of the flow leaving the primary zone 
would predict nitric oxide levels just slightly higher than the 
rapid-quench model represented by the dashed lines. On the other 
hand, the gradual dilution of primary zone products represented by a 
low value of 3 allows nitric oxide to continue being formed at a 
significant rate for a longer time, yielding a higher predicted 
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value, Cii) At overall fuel-air ratios in the vicinity of 0,045 the 

primary zone is considerably richer than the value for peak nitric 

% 

oxide formation. As dilution of the primary zone flow occurs, the 
nitric oxide formation rate at first increases strongly and then 
decreases sharply once diluted below the fuel-air ratio of peak 
nitric oxide formation. In this region, a high 8 value represents 
rapid dilution past the nitric oxide formation peak, while a low g 
value represents a slow dilution toward the peak. Intermediate 
values of 8 allow the dilution process to dwell near the nitric 
oxide formation peak for significant periods of time, yielding high 
nitric oxide predictions. 

A better description of the time history of the mixing intensity 
8 is to realize that it is a function of the intensity of turbulent 
fluctuations within the combustor. Turbulence is generated at the 
entry plane by the flow passing through the blockage plates and 
swirlers and by the combustion in the primary zone. But once the 
flow leaves the primary zone and dilution begins, the turbulence 
decays. One possible assumption for the form of the corresponding 
decay of 8 is 


8 - B exp (-t/Tj (25) 

o a 

where 8 Is the initial 8* t is the time elapsed since exit from the 
o 

primary zone, and is a characteristic decay time, A 8 Q of about 
5000 sec"^ approximates the initial mixing intensity as can be seen 
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from the low fuel-air ratio regime of Figure 11 where the initial rate 

of quenching is the most important factor. The value of 8 will be 

o 

examined later. Figure 12 shows predicted nitric oxide levels for a 
range of values of the characteristic decay time. For comparison, 

the total combustor residence time is 5-6 msec. Note that for a decay 
time longer than this residence time, the predictions are nearly the 
same as for a constant 8 of 5000 as shown in Figure 11. A short 
decay time is qualitatively like the smallest value of constant 8 
case shown in Figure 11. Intermediate values of decay time yield 
results which match the data closely; a decay time of 0.35 msec will 
be used for further calculations. This value of T d suggests that the 
turbulent energy, generated by the swirlers and blockage plates at 
the combustor entrance plane, is dissipated within a distance of 
order one blockage plate diameter. 

From equation (22) it can be derived that for a given combustor 
residence time, the total number of mixing interactions is proportional 
to the time integral of the mixing intensity 8. If equation (25) 
holds, the time integral of 8 is approximately r^8 Q » assuming that 
the residence time is much longer than the decay time For the 

case of rapidly decaying 8 the nitric oxide predictions should not be 
too sensitive to changes in 8 Q and provided the product T d 8 Q is 
fixed. This effect is illustrated in Figure 13 for a T d 8 Q product 
of 1.75. The differences are most marked at overall fuel-air ratios 
around 0.02, where the initial dilution rate is most important. 

Choosing the mixing intensity history 



- 40 - 



fuel - air ratio 


Figure 13. Comparison of N0 X Predictions with 
NASA Swirl-Can Data (4) for Various 
Initial Mixing Intensity/Decay Time 
Pairs (solid curves); fixed primary 
zone and pair product. 



Figure 14. Comparison of N0 X Predictions with 
NASA Swirl-Can Data (4) using 
Parameter Values Fit at Different 
Inlet Temperatures (solid curve) 
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8 * 5000 exp(-t/0.35 msec) sec ^ C26) 

as a reasonable match to emission data at 900 °F inlet temperature, 

the predictions were recalculated for an inlet temperature of 600 °F. 

The results are presented in Figure 14 • The predicted nitric oxide 

levels are close to the actual data, falling somewhat lower near the 

overall stoichiometric fuel-air ratio of 0.069. However, some of the 

model parameters, such as the primary zone mixing parameter s , and 

P 

functional form of B(t), were evaluated only approximately. The 
curve shown in Figure 14 is not necessarily the best match to the 
experimental data. 

Possible mixing model refinements which were evaluated were as 
follows: (i) Calculate the primary zone nitric oxide formation by 

using equation (17) rather than assuming uniform residence times for 
all the fluid elements in the distributions. (ii) Allow the mixing 
intensity to decay toward some non-zero value characteristic of 
"pipe-flow turbulence" rather than decaying toward zero. Neither of 
these modeling refinements produced significantly better predictions. 

A third refinement which was not evaluated would be to attempt 
to reach a given level of mixing parameter s at the combustor exit. 

This value is about 0.5 for the calculations in Figures 13 and 14 which 
is higher than the values of 0.33-0.42 found to give a better match to 
the data in the stirred reactor model for high fuel-air ratios 
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described earlier in this subsection. 


Table 5 presents some data relevant to the NO^ model calculations 
presented in Figure 12 (the "t^ =0.35 msec" curve) and also in Figure 
13 as curve B. 

TABLE 5 

Data Relevant to a N0 X Prediction Curve 
for the NASA Swirl-Can Combustor 


Overall 

Fuel-Air 

Ratio 

Primary 

Zone 

Fuel-Air 

Ratio 

Primary 

Zone 

Residence 

Time 

msec 

Total 

Residence 

Time 

msec 

Number of 
Mixing 
Interactions 

Nitric Oxide 
Mass fraction 

0.0i:> 

0.0375 

0.88 

3.99 

441 

5.69 x 

10“ 5 

0.025 

0.0625 

0.73 

3.24 

440 

1.76 x 

io -4 

0.035 

0.0875 

0.71 

2.88 

439 

3.09 x 

10 _4 

0.045 

0.1125 

0.73 

2.70 

438 

4.53 x 

10~ 4 


inlet temperature = 760 K (907 °F), pressure = 5.8 atm (85.6 psia) , 
mixing intensity as in equation (26). 
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4. Operating Characteristics of the NASA Combustor Inferred from NO Data 

x 

The matching of model predictions to measured NO^ emissions data 
has identified several operating characteristics of the NASA swirl can 
combustor. These can be compared with the requirements outlined in 
Subsection 2.2 for substantially reducing N0 x emission levels below 
current technology values, viz. a fuel lean primary zone, short 
primary zone residence time, significantly more uniform fuel-air mixing 
within the primary zone, and rapid dilution of the primary zone burst 
gases. 

This analysis indicates that the average fuel-air ratio in the 
module wake, the primary combustion zone, is fuel lean for overall 
fuel-air ratios lower than about 0.026. A ratio of about 2.5 between 
primary zone and overall fuel fraction gives a good fit to the data. 

Thus at take-off conditions for current large engines (an overall fuel- 
air ratio of about 0.023) the primary zone equivalence ratio is about 
0.87 which is close to the peak NO formation rate in Figures 4 and 5 
for nonuniform mixtures. Greater entrainment of air into the primary 
zone than occurs in the current combustor design would therefore be 
desirable. 

The values of s found necessary to match the data (about 0.4) are 

also not significantly better than is achieved with current technology 

(3) 

combustors. A previous review v 7 indicated that s values between 0.3 

and 0.7 matched N0^ emissions data obtained from conventional combustors. 
If the primary zone in the NASA swirl can be leaned out further, as 
suggested in the previous paragraph, then Figures 3 through 6 indicate 
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that significantly better mixing in the module wake would improve the N0 x 
emission characteristics. 

The primary zone residence time is 0.7 - 0.9 msec, which is below 

conventional combustor values. Thus, some reduction in NO emissions 

x 

has been achieved through the use of many modules each with its own 
primary zone. 

A comparison of the dashed lines in Figures 13 and 14 with the dilution 
zone model predictions (solid lines) indicates that at overall fuel-air 
ratios of about 0.025 a significant part of the exhaust NO has been formed 
in the dilution zone, downstream of the primary zone exit. The mixing 
model studies suggest that within about one blockage plate diameter of the 
primary zone exit this dilution occurs rapidly. However, the mixing in- 
tensity then decays. Increasing the mixing intensity between dilution air 
and primary zone exit gases would further reduce N0^ emissions. 



Figure 15. N0 X Predictions versus Pressure for a NASA Swirl- 
Can Module, Baseline Case; see text for details; 
circles are test rig data obtained by NASA (16) 
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5. Parametric Studies 

Two parametric studies are presented in this section. The goal of the 
first study is the extrapolation of NC> x emission index data obtained from 
a combustor test rig. In particular, it was desired to extrapolate from 
test rig pressures to the higher pressures typical of engine operating con- 
ditions. For this study the externally imposed factors were as follows: 
air inlet temperature = 733 K (860 °F) , overall fuel-air ratio = 0.02 
(for kerosene fuel), reference velocity = 30 m/s (100 ft/s), pressure 
ranging from 6 to 25 atm. 

Modelling parameters determined in Subsection 3.2 were used as baseline 

values (s =0.4, F/F =2.5, V =8.7 in^, 8 = 5000 sec \ and T = 0.35 

p p o p o a 

msec). Figure 15 shows the baseline predictions plotted against pressure, 
assuming that the modelling parameters are not functions of pressure. The cir- 
cles in this figure are test rig measurements obtained by NASA. The predictions 
are in close agreement with the data in both absolute level and pressure 
trend over the range of available data. This gives some credibility to 
the predictions at the high pressures. The N0 x predictions do not lie along 
a straight line in this log-log plot against pressure. The curve is almost 
linear at low pressures, but the slope decreases at higher pressures due to 
the increasing importance of NO-consuming reactions as equilibrium is 

approached. For the conditions pertaining to this baseline condition, the 

0.59 

pressure dependency of the emission index can be expressed as p * between 

6 and 10 atm, decreasing to p^’^ between 15 and 25 atm. These are com- 

1/2 

parable to the simple p dependency often quoted in the literature, e.g.. 
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1/2 

Heywood and Mikus (3) . The p dependency is based upon a constant 
burned gas temperature, rather than the constant inlet air temperature 
used in this section. 

NASA has experimented with many different swirl-can geometries (16) ; 
NO^ emissions pressure dependency has ranged from p^’^ to p^*^. One 
objective of this parametric study was to see if modest variations of 
the parameter values would significantly change the pressure dependency 
of the NO^ predictions. 

In order to test the sensitivity of the N0 x predictions to vari- 
ations in the modelling parameters for this type of extrapolation, each 
of the five parameters was varied one at a time. The resulting predictions 
are shown in Figure 16. Curves have been drawn for the case of modelling 
parameters which do not change with pressure. The effect of parameters 
which do change with pressure will be discussed later. 

Figure 16a shows that the predictions are rather insensitive to 
small variation in the primary zone mixing parameter s^ for this particu- 
lar baseline. This is due to the mean equivalence ratio (0.75) of the 
primary zone which results from values of the overall fuel-air ratio (0.02) 
and the primary zone enrichment parameter (2.5). As can be seen in 
Figures 3-6, in the vicinity of $ = 0.75, the NO^ formation rate is reason- 
ably insensitive to s^. Note that there would be a great sensitivity at 

$ =0.5, for example, which would be the case for combustor operation at 

P 

an overall fuel-air ratio of 0.013. 

Figure 16b shows the sensitivity of modelling prediction to variation 
in primary zone enrichment factor. For the case being calculated, the 
values F p/ F 0 = 2.75, 2.5, and 2.25 correspond to primary zone equiva- 
lence ratios If! = 0.83, 0.75, and 0.67, respectively. Referring 
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again to Figures 3 through 6, for s =0.4 the primary zone NO formation 

p x 

rate is a very weak function of cj^ . However, in the dilution zone the 
leaner mixtures will be more readily "quenched 11 to negligibly small NO^ 
formation rates. This explains the sensitivity to primary zone enrich- 
ment factor apparent in Figure 16b. 

Sensitivity of the N0 x predictions to primary zone volume is pre- 
sented in Figure 16c. As noted in Subsection 3.2, is proportional 
to the primary zone residence time, and primary zone NO^ is nearly linear 
with this time. However, for the case under consideration, significant 

NO is formed in the dilution zone; independent of V . Thus, the 
x p 

overall dependency of N0 x on is less than linear. This dependency 
is slightly weaker at high pressures, since NO comes closer to equilibrium 
levels. 

Figures 16d and e demonstrate the effects of varying dilution zone 
parameters. For a fixed mixing decay time - 0.35 msec, higher values 
of the initial mixing intensity represent earlier, more complete mixing. 
Since the current case has a lean primary zone, this mixing simply leans 
out the primary zone products, thus quenching N0 x formation. This in- 
fluence of 3^ is shown quantitatively in Figure 16d. The parameter 
has a similar but less marked effect. For a fixed = 5000 sec \ higher 

values of t, represent more complete mixing. However, for a given degree 
d 

of mixing, the mixing does not occur as early as it did for higher va- 
lues of Figure 16e presentes the NO^ predictions for a range of 

v 

None of the curves in Figure 16 exhibit a N0 x prediction pressure 
dependency which is noticeably different from the baseline case. There- 
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fore, experimentally observed pressure dependencies which differ signifi- 
cantly from the trend of the baseline case must be explained in terms of 
pressure-dependent modelling parameters. 


The effects of modelling parameters which vary with pressure can 

be assessed as in the following example. In Figure 16c, imagine a line 

3 

drawn from low pressure end of the V = 10.4 in curve, to the high 

3 

pressure end of the = 7.2 in curve. This would correspond to a 
primary zone whose volume decreases by 31 percent as the operating pres- 
sure increases by a factor of 5. The slope of that imaginary line would 
be 0.41, implying a pressure dependency of p^'^ for the NO^ emissions. 

In this way, a range of pressure dependency exponents can be explained 
in terms of modelling parameters which vary with pressure. Current 
understanding of the five modelling parameters is insufficient to allow 
calculation of how they might change with combustor operating conditions. 


The second parametric study examined operation of the NASA swirl- 
can combustor at conditions typical of state-of-the-art jet engines. 

The modelling parameters were investigated with the goal of improving 
the predicted NO^ emission index. In this case the externally imposed 
factors were: air inlet temperature = 775 K (935 °F) , pressure = 23 atm, 

overall fuel-air ratio = 0.023 (kerosene fuel again), and reference 
velocity = 18 m/s (60 ft/s). 


Once again, modelling parameter values from Subsection 3.2 were 
used as a baseline. The baseline emission index prediction is 38 g 
N 02 /kg fuel. This is comparable to measured values for conventional 
combustors. The emission index should be decreased to approximately 
10 g/kg fuel in order to meet proposed emission standards. It 
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was not possible to predict an emission index less than 10 by reasonable 
variations of a single modeling parameter from the baseline case. 

The modeling parameters were therefore varied in pairs in order to 
explore and take advantage of their interactions. Results are presented 
in Figures 17 through 20. For reference^ the baseline primary zone 

equivalence ratio 4> =0.87. 

P 

Figure 17 presents emission index predictions for ranges of and 
As can be seen in Figure 5, the primary zone NO^ should decrease 
with increasing s^; and due to the lean primary zone, more intense mixing 
will simply quench the N0 x reactions faster. Similar results for ranges 
of Sp and are shown in Figure 18. However, neither of these pairs 
of varied parameters predicted an emission index significantly below 
20. At this mean primary zone equivalence ratio, simply changing 


mixing rates is not the answer to significantly lower N0^ emissions. 

The effects of V and F / F for the current case are presented in 
p p o 

Figure 19. Primary zone enrichment factors F^/F q = 1-5, 2.0, 2.5, 

and 3.0 correspond to = 0.51, 0.69, 0.87, and 1.05, respectively. 

Again referring to Figure 5, for s = 0.4 the primary zone NO formation 

P x 

rates are roughly equal for ^ = 0.87 and 0.69, somewhat lower for 
= 1.05, and approximately a factor of 2 lower for = 0.51. This 
relative ranking is observed in Figure 19 at the high end of the V 


scale, where long primary zone residence times overshadow dilution zone 


processes. For low values of the dilution zone has a higher relative 
significance. The rich primary zone corresponding to F^/F^ = 3.0 con- 
tinues to form N0^ well into the dilution process, thus overtaking the 
leaner primary zone predictions. Varying only and F^/F^, emission 
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indices lower than 10 are achieved only for very lean, small primary 
zones. Unfortunately, this combination is likely to produce an 
unacceptable level of CO due to incomplete combustion (not included 
in this modeling effort). Thus, variations which do not consider mixing 
are also not the answer to low NO in this case. 


More promising predictions are presented in Figure 20, where s^ 
and Fp/F Q are varied. These results are fundamentally similar to Figure 
5, but include the effects of finite residence time and a dilution pro- 
cess. An emission index lower than 10 is predicted in this case for a 

lean, well-mixed primary zone, e.g. F /F £ 1.75 (<p 0.60) for s=0.2. 

p o p 



2 4 6 0 10 20 


Figure 19. N0 X Predictions versus Primary Zone 
Volume for Various Values of Primary 
Zone Enrichment Factor, Typical 
Engine Conditions; see text for 
details 


I 
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6. Air Flow Streamline Calculations 

As a first cut in understanding the flow field downstream of the 
NASA swirl-can modules, an existing finite-difference computer program 
was adapted. This program was TEACH-T developed by Spalding et al, (17), 
which calculates steady, non-swirling, recirculating, axisymmetric , tur- 
bulent flow fields. The boundary conditions were set up as diagrammed 
in Figure 21: a cylinder with cross-sectional area equal to the liner 

area per module at the blockage plate plane of the full NASA combustor; 
zero velocity at the walls; at the blockage plate plane no flow at 
swirler hub radii, uniform axial flow at swirler vane radii, no flow 
at blockage plate radii and uniform axial flow at radii between the 
blockage plate and the cylinder wall; zero pressure gradient 4.5 dia- 
meters downstream of the blockage plate; and 19 present of the mass flow 
passing through the swirler (see Subsection 9.2). Within the flow 
field gas properties were taken to be those of ambient air. 

Figure 22 is a streamline plot derived from the computer calculation 
for 38 m/s (124 ft/s) mean axial velocity. The flow is from left to 
right. The streamlines toward the top of this figure represent the 
annular jet issuing from the gap between the blockage plate and the 
cylinder wall. Four streamlines can be seen emerging from the 
swirler; they move rapidly outward, so that the angular velocity of 
this flow in the real system would be lost quite early. The swirler 
streamline closest to the axis of symmetry (bottom of the figure) is 
the contour of zero stream function, and thus outlines the mean recir- 
culation zone. This zone extends about 1.8 blockage plate diameters down- 
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Figure 20. N0 X Predictions versus Primary Zone 
Nonuniformity for Various Values of 
Primary Zone Enrichment Factor, 
Typical Engine Conditions; see text 
for details. 



dimensions in inches 


Figure 21. Schematic of Boundary Conditions 
for Numerical Flow Calculations. 
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stream and about 13% of the total flow is involved in recirculation. 

It must be noted that these are the mean streamlines calculated for 
this turbulent flow field; they do not show the turbulent mass transfer 
which crosses such streamlines. 

Streamline calculations were also performed for 9 m/s (31 ft/s) 
axial velocity, with virtually no change in the streamline pattern. 

This represents a range of Reynolds numbers of 35,000 to 140,000, based 
on blockage plate diameter. From these calculations it follows that the 
size of the module wake primary zone, and the fraction of the air which 
passes through it are weak functions of combustor operating conditions. 
This supports the choice of and F^/F^ as primary zone parameters 
which are fixed for a given combustor geometry. In addition, the 

3 

value chosen for (8.7 in ) is of the right order. 



Figure 22. Streamline Plot of Calculated Flow 
Field; flow is from top to bottom, 
with symmetry axis at left side 
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7. Calculation of Fuel Droplet Evaporation 

As argued in Subsection 2.2, the air flow pattern and early fuel 
processes are of fundamental importance for determining NO^ formation. 
Section 6 dealt with the calculation of air flow; this section will 
address the fuel processes. 

As will be seen in Subsection 8.2, fuel droplets are formed by 
the breakup of sheets and ligaments associated with the swirler vanes 
The maximum stable droplet size may be calculated as suggested by 
Prandtl (18, p. 326) 


d % 15.4 — (27) 

pv 2 

where a is the fuel T s surface tension, p is the air density, and v 
is the relative velocity between fuel and air. This equation will be 
verified in Subsection 8.2. For a given fuel a is a function of 
temperature, so it is necessary to know the temperature of the fuel 
droplet; this was assumed to be the wet bulb temperature, or 0.99 
times the absolute boiling temperature, whichever was lower. The 
quantity pv 2 was calculated for inlet conditions using the open area 
at the plane of the blockage plate. 

A computer program was written to calculate the evaporative 
flight of a fuel droplet. Komiyama (19, Appendix 3) described a 
similar program. Stagnant droplet evaporation rates based upon 
Goldsmith and Penner (20) were corrected for forced convection 
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following Ranz and Marshall (21). Droplet acceleration from an 
Initial state of rest was calculated using the standard drag equation 
for a sphere, corrected for evaporation effects according to 
Spalding (22). Calculations were conducted for projected idle and 
maximum cruise power conditions; evaporation would be somewhat faster 
at take-off. Conditions for the gas surrounding the fuel droplets 
were taken to be alternately inlet conditions and combustor exhaust 
conditions. Even the latter assumption is rather conservative since 
it is far below the combustion temperatures expected in the flame 
zone. Table 6(a) presents the results of these calculations. The 
higher gas temperature assumptions predict essentially complete 
evaporation in less than 0.1 module diameter, even at idle. This 
supports the assumption of instantaneous evaporation which was made 
in the NO^ model. Note, however, that at idle a fuel droplet would 
evaporate relatively slowly in regions not heated by combustion 
products. 

Similar calculations were done for an ambient pressure and 
temperature simulation which will be described in subsection 10.2. 

These results are shown in Table 6(b). It was hoped that pentane, 
which has a normal boiling point of 309 K, would be sufficiently 
volatile to simulate the evaporation of kerosene at the higher 
temperatures to which it would be exposed in an engine’s combustor. 
However, due to lower gas density and somewhat higher surface tension, 
the pentane droplets are significantly larger [see equation (27)] . 

For a nominal gas temperature of 1000 K, the pentane requires more than 
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one module diameter to evaporate, even at simulated cruise conditions. 
This explains the experimental difficulties which will be discussed 
in subsection 10.2. 


TABLE 6 

Fuel Droplet Evaporation Parameters and Calculated Results for Engine 

Conditions and a Simulation 


(a) kerosene fuel 
estimated 

idle cruise 


pressure, atm 

4. 

2 

13, 

2 

1 

1 

gas temperature, K 

482 

1005 

763 

1533 

1000 

1000 

fuel temperature, K 

441 

524 

580 

617 

306 

306 

gas velocity, m/sec 

49 

49 

62 

62 

49 

64 

initial drop size, y 

21 

13 

22 

11 

75 

43 

99% evap. time, msec 

5.95 

0.16 

0.56 

0.04 

5.3 

1.9 

distance, module d's 

4.99 

0.07 

0.52 

0.01 

2.3 

1.1 


(b) pentane fuel 
simulated 
idle cruise 
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8. Flow Visualization Experiments 
8.1 Flow hardware 

A program of flow visualization experiments was undertaken in 
order to obtain a better qualitative understanding of the flow field 
in the NASA swirl-can combustor. A stainless steel blockage plate 
was made with area equal to the average area of the 120 blockage 
plates in the full NASA combustor. This plate was silver soldered to 
a standard swirler. A plexiglass can body was fashioned and pressed 
onto the swirler. This assembly was mounted by two struts, centered 
in and flush with the downstream face of a thick plexiglass flange 
whose open area was equal to the liner area per module in the full 
NASA combustor. The flange was bolted to the downstream end of a 
slightly larger plexiglass tube. A fuel tube was fashioned such 
that it terminated 1/8 inch upstream of the swirler with a 3/16 inch 
tube having an 0.052 inch orifice at the end; this tube could be 
aimed at any part of the swirler during the experiments. Details 
of the visualization test section are shown in Figure 23. 

Air was supplied to this test section from a 200 cfm compressor 
through a pressure regulator, throttling valve, ASME orifice rig, 
flow straight ener, and 20 diameters of straight pipe. The test 
section generally discharged directly into the room, and ambient 
condition prevailed at the module. The air flow rate was controlled 
so that reference velocities could match those of the full NASA 
combustor. Reynolds numbers based on blockage plate diameter were 
of order 10^. 



dimensions in inches 


Figure 23. Diagram of Visualization Test Section 
flow is from top to bottom 



Figure 24. Helium-filled Soap Bubble Traces in 
Module Wake; mean flow is from top 
to bottom; bubbles are injected at 
tip of recirculation zone 
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8.2 Techniques and Results 

The size of the recirculating wake downstream of the module 
was determined by probing with a wool tuft. Probing along the axis 
of the flow field, a length of 1.0 module diameter was observed; 
this increased to 1.8 diameters when a confining tube section was 
added downstream of the module; no Reynolds number effect was observed. 
These results quantitatively confirm results of the flow calculations 
presented in Section 6; they also support the value chosen for the 
primary zone volume in 9ubsection 3.2, as well as the assumption that 
Vp is a weak function of operating conditions. The remainder of the 
experiments described in this subsection were conducted without 
the downstream confining tube mentioned above, to allow easy 
photographic access. The size of the module wake region was affected 
by this omission, as described above, but the flow was expected to 
remain qualitatively the same. 

A smoke generator was connected to the fuel tube to visualize 
the flow of a gaseous fuel. Running the test section at less than 
about one tenth nominal air flow, the smoke density was sufficient 
to see a clearly defined module wake "bubble” with smoke distributed 
rather evenly throughout. Downstream of this bubble the smoke was 
considerably less dense. These results support the assumption of a 
s tirred-reactor primary zone which operates considerably richer than 
the overall fuel-air ratio, as proposed in Subsection 2.2. 

An alternative technique for visualizing streaklines in the flow 
is the addition of helium-filled soap bubbles. A Sage Action Inc. 



62 


Model 3 bubble generator was used in a setup which paralleled the 
smoke visualization. The size of these bubbles could be adjusted 
such that they were neutrally buoyant in air, thereby following the 
air flow patterns with no inertial effects; this size was approximately 
0.1 inch, less than 0.05 module diameter. Unfortunately, these 
bubbles burst rather than flow through the swirler vanes, due to the 
narrow passage and high shear. Therefore, bubbles were introduced 
into the module wake region either along the outer edge or, as shown 
in Figure 24, at the downstream tip along the axis. The mean air 
flow in this figure is from top to bottom at about one fifth nominal 
flow rate. Exposure time was 1/50 second, about one residence time, 
with forward scatter lighting such that a bubble is seen as two 
flares separated by one bubble diameter. In this picture, the 
bubbles are seen to travel "upstream" toward the swirler, then 
undergo looping motions indicative of the recirculation pattern, 
becoming involved in tight eddying motions along the way. To follow 
the bubble motions, high speed motion pictures at 5,000 frames per 
second were also taken. Projected at 16 frames per second, these 
movies show a module wake region full of bubbles "lazily" wandering 
around. There does not appear to be any organized flow pattern except 
for a tendency to move toward the swirler when a bubble is near the 
axis, and away from the swirler when near the outer edge of the 
module wake region. One is reminded of a random-walk process. There 
is never a coherent motion of a significant number of bubbles, or any 
evidence of large eddy shedding. When viewing the experiment 
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directly, the eye integrates these motions into a vibrating toroidal 
vortex. This helium-filled soap bubble experiment gives additional 
support to treating the primary zone as a stirred reactor. 

To visualize the fuel droplet processes in the extreme limit of 
a nonevaporating fuel, a set of experiments was conducted in which 
water was introduced through the fuel tube. Figure 25 is a 1/20 
second exposure at about one half nominal air flow rate, with a 
water-air ratio of 0.027 (comparable to the take-off fuel-air ratio 
of an engine). The air flow is from right to left. One can see the 
water drops coming out of the swirler and flowing outward as they move 
downstream; the details are out of focus except near the bottom of the 
picture. These droplets appear to be following an air flow pattern 
similar to the streamlines predicted in Section 6. Note, however, 
that the droplets do not turn back inward as those streamlines 
would predict. Apparently, the inertia of the water droplets carries 
them across the air streamlines into the annular air jet which comes 
around the blockage plate. This jet shatters the droplets down to 
a smaller size which is difficult to follow photographically. These 
droplet processes were confirmed by stroboscopic pictures using a 
3 microsecond flash, and also with motion pictures taken at 5,000 
frames per second. The movies enable one to follow individual 
droplets from frame to frame as they go through the processes described 
above. The process by which the droplets are formed coming out of the 
swirler was studied by stroboscopic close-up pictures such as 
Figure 26, taken at nominal air flow and a water-air ratio of 0.035. 
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Figure 25. Streaklines of Water Droplets in 
Module Wake; flow is from right 
to left; water is injected behind 
swirler . 



Figure 26. Closeup of Water Emerging from 

Swirler; water is injected behind 
swirler at center of hub. 
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Here can be seen sheets of water on the swirler vanes, ligaments 
being tom off, and droplets in various stages of development. 

All of the water experiments thus far described were photographed 
using forward scatter lighting. It is difficult to obtain quantitative 
information about droplet sizes from this method because of the lens-like 
behavior of the water drops. So a series of shadowgraph pictures 
were taken by lining up the strobe light, a Fresnel lens, the field 
of droplets, and the camera along the same optical axis. Similar 
techniques have been described by Frazer and Dombrowski (23). 

Figure 27 is a shadowgraph picture taken parallel to the blockage 
plate at nominal air flow with a water-air ratio of 0.039. Air 
flow is from top to bottom, and only the left side of the flow field 
is shown. Measurements from photos such as this one demonstrate the 
effect of reference velocity on droplet size. Data are presented in 
Table 7, along with predictions following Prandtl, equation (27), 
and Wigg (24). The Prandtl equation apparently predicts an 
intermediate droplet size, and was therefore used to estimate initial 
droplet sizes for the evaporation calculations presented in section 7. 

The Wigg correlation is in poor agreement with measured droplet sizes, 
having been developed for more conventional fuel injectors. 

A final set of visualization experiments were carried out under 
burning conditions. For these, the downstream face of the plexiglass 
flange in which the module was mounted needed to be insulated from 
the combustion. A 1/4 inch thick transite sheet was fashioned with 
a hole to match the plexiglass flange hole. The assembly was 
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Figure 27. Shadowgraph of Water Droplets in 
Module Wake; flow is from top to 
bottom; only left side of flow 
field is shown 


TABLE 7 

Measured and Predicted Droplet Sizes 


air reference velocity, m/sec 13.5 19.2 

representative drop measurements, y 

near swirler 850 470 

far downstream 240 120 

predicted drop sizes, y 

Prandtl, equation (27) 440 219 

Wigg ( 24; 82 58 


38.4 

90 

mist 
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mounted on top of an existing vertical chimney more than one foot 
square in internal dimensions. Burning experiments were conducted 
using propane (gas) and pentane (liquid) fuels, at reference 
velocities and equivalence ratios typical of take-off and idle 
conditions. Figure 28 shows the combustion of propane at an 
equivalence ratio of 0.21 (approximately idle), for a reference 
velocity about 40% of idle conditions. The exposure time was 1/50 
second, using illumination from the flame only. Flow is again from top 
to bottom. A roughly annular luminous region can be seen extending 
from the swirler about one module diameter downstream. It is a blue 
flame. Behind the flame are two light bands with a dark band between 
them, these are the edge of the hole in the transite sheet which 
protects the plexiglass flange, and defines the outer limit of the 
annular air jet. The visible flame was similar for higher reference 
velocities, becoming longer when operating richer. The burning of 
pentane did not visibly differ from the propane results described 
above, except that at low air flow rates a flickering yellow flame 
could be seen near the burner axis, Immediately downstream of the 
swirler. This secondary flame was attributed to the large droplets 
of pentane which can exist at low air velocities [see equation (27)] . 
However, even under simulated take-off conditions, with <(> ^ 0.4 some 
pentane droplets did escape from the combustion zone. These droplets 
were made visible by the two stroboscopic techniques described above, 
shadowgraph and forward scatter lighting. An example of the latter 
method is presented as Figure 29, which shows the right side of a 


Figure 28. Combustion of Propane in Module Wake; 

overall equivalence ratio - 0.21; 
flow is from top to bottom, fuel is 
injected behind swirler 



Figure 29. Combustion of Pentane in Module Wake 
with Droplets Escaping; overall 
equivalence ratio = 0.41; flow is 
from top to bottom, fuel is injected 
behind swirler; droplets illuminated 
by strobe light 
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top to bottom flow field. The lens open time for this picture was 
1/200 second, forming an image of the blue flame at the right side 
where the annulus is tangent to the line of sight, thus producing 
the greatest intensity of light, While the lens was open, a 3 
microsecond flash illuminated the pentane droplets, freezing their 
motion. Although it seems that a significant amount of fuel might 
be escaping combustion, in the full NASA combustor there would be 
an adjacent module wake to assist in the burning of these stray 
droplets; in addition, droplets would vaporize much more readily in 
the full NASA combustor, as demonstrated in Section 7, 

The visualization experiments qualitatively confirmed many of 
the simplifying assumptions about the flow field which were made in 
Subsection 3.2. 



Figure 30. Diagram of Test Section for Cold 
Flow Tracer Experiments; flow is 
from left to right 
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9. Nonburning Tracer Measurements 

9.1 Flow hardware and sample probe 

Quantitative data about the flow in the NASA combustor was 
obtained by taking pressure measurements and cold flow tracer gas 
samples from a single-module test section, shown in Figure 30. The 
flow areas of the axisymmetric test section are an idealization of 
the flow areas existing in the full annular combustor liner, divided 
by the number of modules. Thus, there is a constant diameter section 
upstream of the module, extending about 0.2 module diameter downstream, 
a 15° included angle conical expansion for about one diameter, and 
then a constant diameter to the end of the test section. A blockage 
plate and standard swirler as described in Subsection 8.1 were silver 
soldered to a standard can body, which was supported in the test 

section by three radial struts. 

To achieve more spatial uniformity of tracer upstream of the 

swirler, a tube was attached to the can body extending it about six 
can diameters upstream. Just past the entrance to this extension 
tube a coarse-mesh screen was installed to support the tracer tube 
and increase mixing. The tracer was injected from a point on the 
axis immediately downstream of this screen, discharging through 
four radial holes^and boosted with compressed air. This additional 
air was less than five percent of the estimated flow through the 
extension tube. 

The test section exhausted into a large pipe tee and then down 
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Into the building exhaust trench- This arrangement provided access 
for a straight gas sampling probe to extend into the test section, 
parallel to the axis. Expressed as blockage plate diameters, the 
probe consisted of about 3d of 0.06d tubing, flowing into an 0.2d 
tube which extended beyond the test section. The probe was mounted 
on a three-dimensional traversing mechanism such that any point of 
the downstream flow field could be sampled. 

Propane was used in low concentrations as the tracer gas. The 
gas sample from the probe was passed through a Beckman Model 400 FID 
hydrocarbon analyzer. Sample flow rate could be varied and was found 
to have little effect on the relative tracer concentrations which were 
observed. 

Pressure taps were installed to monitor pressure inside the can 
body and upstream and downstream of the blockage plate. An upstream 
total pressure probe was installed when the tracer tube and can 
extension tube were not in place. 

9.2 Technique and results 

The pressure measurements obtained from the locations described 
above were used to calculate that approximately 19 percent of the air 
flow passes through the can body. This calculation was derived from 
compressible flow relationships and the known areas, assuming 
isentropic flow upstream of the blockage plate and swirler. The open 
area of the swirler is 21 percent of the total open area in the plane 
of the blockage plate. This indicates that the initial velocities 
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of the swirler jets and annular air jet are nearly equal. 

With the extension tube and tracer tube in place, gas samples 
were withdrawn from the flow field at several axial stations downstream 
of the blockage plate. After determining that the tracer concentrations 
were symmetric, the probing was restricted to one quadrant of the 
test section. Probing was done on a "grid point" basis. The grid 
was as fine as 0.050 x 0.100 inch (0.022 x 0.045 blockage plate 
diameter), depending upon the steepness of concentration gradients. 

A three-dimensional representation of the tracer concentrations 
measured at 0.045 diameters downstream of the blockage plate is 
presented in Figure 31. The two axes at the left of this figure 
correspond to the orthogonal radii which bounded the probed quadrant. 

The axis extending to the right corresponds to the axis of the 
test section; distance along this axis represents tracer concentration. 
The flat region near the axis has a concentration equal to the 
overall tracer to air ratio. Further out one sees two complete and 
two partial concentration peaks which represent the rich jets coming 
from the swirler; since there are twelve vanes in the complete 
swirler, three jets can be expected in each quadrant. Maximum 
concentrations are about seven times the overall level, roughly 
corresponding to the flow split estimate described earlier. Beyond 
the peaks can be seen an outer plateau where concentrations are 
between three and four times the overall value; this represents the 
relatively stagnant region immediately downstream of the blockage 
plate (see the streamlines of Figure 22), where fuel and air have 
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been mix ed and a flame may be stabilized. The combustion observations 
described in Subsection 8.2 support this location of a stabilization 
region. 

A more quantitative presentation of the tracer concentrations can 
be had from the contour maps of Figure 32. These maps correspond to 
axial stations at 0.0045, 0.045, 0.09, 0.18, and 0.45 diameters 
downstream of the blockage plate. Each map represents one quadrant 
of the flow field, viewed from downstream. The test section axis 
is at the upper left, with orthogonal radii running along the top 
and left edges. The arc running from upper right to lower left 
lies slightly inside the test section wall; tracer concentrations 
outside this arc are essentially zero. Each boundary between a 
shaded and an unshaded region is a contour line (constant tracer 
concentration) . Within each map the contour interval is one tenth 
of the maximum conentration at that axial station. These maxima are 
7.0, 6.2, 4.9, 3.1, and 1.5 times the overall for the five maps. In 
each map, the region nearest the axis (upper left) is at the overall 
concentration; concentrations increase with radius, passing through 
the rich jet region, then decreasing toward zero as the test section 
wall is approached. The jet structure is clearly defined in the 
first few maps, becoming more diffuse and merging circumferentially 
as the flow progresses downstream, A small amount of counter- 
clockwise swirl is also apparent. By 0.9 diameter downstream of the 
blockage plate, the tracer concentration is quite uniform across the 
test section. Decreasing the reference velocity by a factor of four 


74 



figure 31- 


1 Plot of Tracer 
Thr ee Dimensions Axial 

Concentrations a defi ne probed 

Station*, axes along horizontal 

distance aJ - ons ._„^ rr atioi 




Contour Map of Tr ^ c : e ^ ^ C Downstream; 

F «"'“ 3 at 0.0045 «° dul "- s “i s a t «PP« left > 

S&STS- ^ — a ” 


75 



Figure 32b. Contour Map of Tracer Concentrations 
at 0.045 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 



Figure 32c.. Contour Map of Tracer Concentration 
at 0.09 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 
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Figure 32d. Contour Map of Tracer Concentrations 
at 0.18 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 



Contour Map of Tracer Concentrations 
at 0.45 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 


Figure 32e 
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did not noticably alter these contour maps, thus confirming the lack 
of Reynolds number effect predicted in Section 6. 

The swirler jets can be quantitatively followed as they progress 
downstream, moving outward and entraining surrounding gases. A 
radial line drawn through the concentration peak of a jet at each 
axial station produces a cross-section representative of the jet 
structure. A number of these jet cross-sections are shown in 
Figure 33, where tracer concentrations are plotted against radius. 

The jets can clearly be seen to move outward and dilute themselves 
with surrounding gases as they progress downstream. This process is 
completed by about one module diameter downstream. 

A gross dilution rate for the swirler jets can be estimated by 
using this tracer concentration data. Following Hoult and Weil (25) 
an entrainment parameter may be defined as 

a “ p A V(1 - cos 6) 

where m is the mass entrainment rate, p is the gas density, A is the 
area through which mass is being entrained, and V(1 - cos 0) is the 
relative tangential velocity between the swirler jets and the annular 
air flow between the swirler jets and the test section wall. All 
dilution is attributed to tangential entrainment for the purposes of 
this calculation. Let d^ be the blockage plate diameter and d Q be the 
test section diameter at the blockage plate axial station. If the 
swirler jets entrain all of the annular air flow, one may write 
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» - I p v < d o - aj) 


C29) 


If the entrainment occurs in an axial distance equal to the blockage 
plate diameter, the entrainment area may be approximated as 


A ^ 

A 'b 


TT d 2 

R 

cos 0 


(30) 


Substituting equations (29) and (30) into equation (28), one obtains 


Or 


a /v 


(d 2 - d 2 ) cos 0 
9 E 

4d 2 (1 - cos 0) 


(31) 


Using a value of 0 % 60* obtained from Figure 33 and plugging in the 

values of d and d , one obtains a value of a ^ 0.08. This crude 

op 

estimate agrees surprisingly well with the value 0,078 (= 0.11/ Jl) 
reported by Hoult and Weil ( 25 ) for Gaussian velocity profiles. 

Cold flow data presented in this section generally support the 
flow calculations of Section 6. Additions have also been made to the 
qualitative understanding of the flow field, particularly the region 
within one diameter of the blockage plate. A quantitative under- 
standing of the flow field under combustion conditions will be 
developed in the next section. 



TRACER CONCENTRATION, 



igure 33. Tracer Concentration versus Radius at 
Various Axial Stations; radius runs 
through swirler jet concentration 
peak at each station 



Figure 34. Diagram of Test Section and Probe for 
Gas Sampling under Burning Conditions; 
a modification of the cold flow tracer 
section (Figure 30); flow is from 
left to right 
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10. Measurements under Burning Conditions 

10.1 Flow hardware and sample probe 

The measurements taken under burning conditions were done in the 
same test section as the cold-flow tracer experiments, as described 
in Subsection 9.1, with the following changes (see Figure 34): 

The upstream extension tube was removed from the swirl-can. The 
tracer injector was replaced by a fuel tube which simulated NASA's 
design, a 3/16 inch tube running along the axis of the test section 
to a point 1/8 inch upstream of the swirler, terminating in an 
0.052 inch orifice. A single-electrode igniter was added Immediately 
downstream of the blockage plate, at the top of the test section such 
that the spark would arc between the igniter and the blockage plate. 

Copper tubing was wrapped around the outside of the test section 
for cooling water. A spray bar was installed inside the exhaust tee, 
downstream of the test section, for water-cooling the walls of the 
tee. A spray-cooled deflection plate was installed in the building 
exhaust trench below the exhaust tee. 

A water-cooled sample probe was constructed of three concentric 
stainless steel tubes, with a skewed copper tip. Mounted on the 
three-dimensional traversing mechanism, the probe passed through the 
exhaust tee, extending into the test section parallel to the axis. 

The portion of the probe which extended into the test section had 
a diameter of 0.20 inch, which is 0,09 times the blockage plate 
diameter. The probe tip (see Figure 34) was contoured to provide an 
aerodynamic ally smooth, transition from the concentric tubes to the 
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sample inlet * This inlet was located approximately one probe 
diameter beyond the end of the concentric tubes, offset from the 
probe* s axis by one radius. Thus, the probe could be used to sample 
right up to the wall of the test section without excessively 
disturbing the flow field. 

10.2 Choice of propane fuel 

As described in Section 7, it was determined that pentane fuel 
would have a longer evaporation distance for ambient inlet air 
conditions than kerosene has under typical engine operating 
conditions. Gaseous propane fuel would simulate the instantaneous 
evaporation of a liquid fuel. A comparison of the combustion of 
pentane and propane could confirm the argument of section 7 that 
droplet evaporation effects are not significant in the NASA 
swirl-can combustor. 

When the experimental burner was operated on pentane, it was 
found that a significant portion of the fuel would reach the 
relatively cool wall of the test section without evaporating. This 
liquid fuel flowed downstream along the wall to the end of the 
test section, where the hot gases dump into the exhaust tee. At 
this point a secondary flame was attached, burning inside the exhaust 
tee. While this situation was an interesting combustion phenomenon, 
it could not simulate the combustion downstream of the 120-module 
NASA swirl-can array. That is, a solid wall cannot simulate 
imaginary surfaces of flow symmetry under conditions where fuel droplets 
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Impinge upon that wall. It was decided to proceed with combustion 
experiments using gaseous propane fuel to approximate the rapid 
evaporation of kerosene expected under engine operating conditions. 

10.3 Gas sample analysis hardware 

The gas samples drawn through the probe described in Subsection 

10,1 were pulled through a length of unheated plastic tubing into a 

diaphragm pump, then through an ice-bath moisture trap followed 

by desiccant tubes, and into the continuous-reading instrumentation. 

These instruments were: a Bedonan Model 315A NDIR for CO^, Beckman 

Models 315A and 864 NDIR f s for CO, a Scott Model 150 paramagnetic 

for 0^, a Scott Model 215 FID for total hydrocarbons (HC), and a 

Thermo-Electron Model 10A chemiluminescent for NO . For a sample 

x 

with more than 100,000 ppm HC (as C^) * the sample was diluted with 
nitrogen in order to bring it within the range of the HC instrument. 
Since the NO^ readings were always less than 50 ppm, it was treated 
as a minor species whose magnitude did not affect the calculated 
values of local equivalence ratio and mixture nonuniformity which 
were evaluated in the combustion experiments. 

10.4 Reduction of data to local equivalence ratio 

When a gas sample is taken from a point in the test section, the 
gas analysis instrumentation described in the previous subsection 
provides time-average measurements of mole fractions of C 02 > CO and 
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C >2 and mole fraction times average carbon number for hydrocarbons, 
all on a dry basis. This information alone is not sufficient for the 
determination of local equivalence ratio, unless some assumptions are 
made. Many of the published methods for calculating equivalence 
ratio, such as the Spindt method (26), rely on a ff water-gas ,? 
equilibrium constant. This assumes that H^O and CO are in a 
quasi-equilibrium with CO^ and H^. While this assumption may be 
reasonably accurate for piston engine exhaust gases, it is unlikely 
to be valid for gas samples drawn from arbitrary points in a complex 
combustor flow field. An alternative assumption was therefore 
developed. 

For a combustor fueled by a light stable hydrocarbon such as 

propane, it may be assumed that the hydrocarbon to which the 

hydrocarbon analyzer is responding is the same as the fuel which is 

fed into the combustor. The justification for this assumption is that 

in regions where hydrocarbon concentrations are high, the fuel is 

likely to be nearly unburned and unreacted; in regions of low 

hydrocarbon concentration it is a minor species and the assumption 

has little effect on the equivalence ratio. For the hydrocarbon 

fuel C H , the average reaction for gases sampled at a point in the 
n m 

flow field may be written for major species as 

♦c H + C»W4)0 2 + 3.76taW4)H 2 - TO.,, g + * Oj + CO 

n m 2 

♦^COj “2 *=\o H 2° +JC H 2 E 2 + I S 2 V 


( 32 ) 
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where $ is the equivalence ratio, T is the number of moles of 
products and :x^ is mole fraction of species i. Balance equations 
may be written for C, H, 0, and as follows: 


n * * TCnX C H + X C0 + 

x„„ ) 
CO' 

(33) 

n m 

2 


m* = T(mx + 2x 

+ > 

(34) 

n m 2 

2 


2(n + m/4) = T(2x + 

x co + 2x C 0 2 + X H 2 0^ 

(35) 

2 



3*76(n + m/A) = T x N 


(36) 


The four dry basis species measurements may be equated to / (1-x^ q) » 

% C0 /(1_X H„0 ) ’ X 0 o /(1_X H o 0 ) ’ and n X C H /(1_X H o 0 ) - Xt ±S alS ° 

2 2 2 n m 2 

known that the mole fractions on the right-hand side of reaction (32) 

sum to unity. Thus there are nine equations written in terms of 

n, m, 4>, T, and the seven mole fractions, x_^. For a given fuel n and 

m may be determined, and the system of equation may then be solved. 

In particular, the equations may be written in terms of the unknowns 

<J>, T» and T x^, resulting in linear equations. 

The method described above was applied to data obtained from the 

test section described in Subsection 10,1, burning C^Hg. 

resulting equivalence ratios and mole fractions seemed reasonable, 

except that for about one-third of the data x__ was negative by a 

H 2 

few percent (in rare instances x Q was negative) . This was 
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attributed to a relatively minor species being quite sensitive to 

inaccuracies in the gas analysis data. The calculated equivalence 

ratios were much less sensitive. 

In order to compensate for these measurement errors, a method 

was devised to solve the system of equations in a least squares 

sense with respect to the measurements, constraining x and x„ 

H2 0 

to be non-negative. The method was built around an available 
routine for the minimization of a sum of squares. The equations 
were solved for the remaining variables In terms of x__ , x__ 4>, 

*h 2 h 2 o 

and T. The minimization routine searched over these four variables 
at each point a feasible solution (for the x^) was calculated. An 
objective function could then be expressed as 


E 

i 


R i (1 “ x H20,c) 


i»c 


1 


(37) 


where is the species i measurement expressed as a dry basis mole 
fraction (dividing by n in the case of hydrocarbons), and the 
subscript c refers to calculated values. The minimization routine 
directed the search such that expression (37) was minimized. The 
point at which this occurred was taken to be the best fit to 
experimental data. Equivalence ratios resulting from this method 
will be used in this section. 
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10.5 Reduction of data to local nonuniformity parameter 

Burner gas sample concentrations may also be used to generate 
a measure of the mixture nonuniformity at the point where the sample 
is taken. Assume that the sample probe pulls in a rapid succession 
of discrete packets of various fuel fractions, and that each packet 
is completely homogeneous. Further assume that each packet is in the 
chemical equilibrium corresponding to its fuel fraction, the inlet 
air temperature, and the fuel properties, i.e. combustion is complete 
with no net heat transfer. Flammability limits may be imposed 
upon this assumption; for the calculations presented herein they 
were arbitrarily set at <J>=0 and <J>=3, wide enough not to significantly 
affect the results. If the fuel fractions of these packets are 
distributed according to equation (15) t average mole fractions of 
CO 2 , CO, 0^, and 1^0 for the mixture may be obtained by calculating 
equilibrium moles of each species per unit mass, for a range of fuel 
fractions, then weighting these on a mass basis using equation (15). 
Average mole fractions of CO^, CO, and 0^ on a dry basis may in turn 
be computed; these are what would be measured on the gas analysis 
setup described in subsection 10.3 where the packets would be quenched, 
mixed, and desiccated on their way through the flow system. 

The only unknowns in the calculations described above are the 
variables F and a in equation (15) , Given a set of gas sample 
measurements, F is best computed from the equivalence ratio method 
of the previous subsection. The sum of squares minimization routine 
was then used to search over the variable 0; for each value, a 
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feasible set of dry-basis CO^, CO, and 0^ mole fractions was 
generated as described above. An objective function was then 
calculated as 


E 

i 



(38) 


using the notation of expression (37), Hydrocarbons were not 
included in expression (38) because they would be sensitive to the 
choice of realistic flammability limits. The minimization routine 
then directed the search over a such that expression (38) was 
minimized. The point at which this occurred was taken to be the 
best fit to the experimental data, given all the assumptions made 
earlier in this subsection. Values of the mixing parameter s could 
then be calculated using equation (16); these will be presented in 
the next subsection. 

The assumptions upon which this subsection is based are 
equivalent to the corresponding assumptions which were made in 
Subsections 2.1 and 2.2 for the calculation of NC> x ; the s parameter 
values which are thereby calculated are the correct ones to 
compare with s values which enter into the NO^ predictions. To the 
extent that the assumptions are in error, the errors will self- 
compensate. 
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10,6 Results 

The hardware described in Subsections 10,1 and 10,3 was used to 
obtain local species concentrations while propane was being burned 
in the test section. After determining that the flow field was 
roughly symmetric in terms of these concentrations, detailed probing 
was conducted in a single quadrant. A grid-point system as fine as 
0.100 x 0.100 inch CO. 045 x 0.045 blockage plate diameter) was used, 
depending upon the steepness of concentration gradients. Full 
quadrants were probed at 0.22, 0.45, and 0.90 diameter downstream of 
the blockage plate, for a reference velocity of 23 m/s (75 ft/s), at an 
overall equivalence ratio of approximately 0.4. Concentration data 
thus obtained were reduced to grid-point values of equivalence ratio 
and mixing parameter s, using the methods described in Subsections 
10.4 and 10.5. 

Equivalence ratio contour maps (similar to the maps of Figure 32 
described in Subsection 9.2) are presented in Figure 35. The 
contour interval for all of these maps is 0.2 (j) unit. The unshaded 
region in the upper left of each of these maps (closest to the 
test section axis) is the contour band 0.6 < <f> < 0.8, or 1.5 to 
2.0 times the overall <j>, In Figure 35a (j> increases with radius to 
a mATciTimm of 2.1 in a region where some rich jet structure is visible; 
beyond this, <J) decreases toward a value of zero near the test 
section wall (the annular air jet), Figure 35b shows <J) increasing 
with a radius to a region of 0,8 < <f> < 1.0, with a pocket of <j) 
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Figure 35a. Contour Map of Equivalence Ratio at 
0.22 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 



Figure 35b. Contour Map of Equivalence Ratio at 
0.45 Module Diameter Downstream; 
test section axis is at upper left, 
quadrant viewed looking upstream 
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Figure 35c, Contour Map of Equivalence Ratio at 0.09 
Module Diameter Downstream; test section 
axis is at upper left, quadrant viewed 
looking upstream 



Figure 36a. Contour Map of Nonuniformity Parameter at 
0.02 Module Diameter Downs tream; test sec- 
tion axis is at upper left, quadrant viewed 
looking upstream; values not meaningful in 
region of mottled shading 
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Figure 36b. Contour Map of Nonuniformity Parameter at 
0.45 Module Diameter Downstream; test sec- 
tion axis is at upper left, quadrant 
viewed looking upstream; values not meaning- 
ful in region of mottled shading 



Figure 36c. Contour Map of Nonuniformity Parameter at 
0.9 Module Diameter Downstream; test sec- 
tion axis is at upper left, quadrant 
viewed looking upstream 


I 
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slightly greater than 1.0 at the upper right; <p then decreases toward 
zero at the test section wall. In Figure 35c, except for pockets of 
<J> slightly greater than 0.8 at the upper right, <f> simply decreases 
toward zero beyond the central region. 

There are significant differences between these (p maps obtained 
under combustion conditions and the cold-flow tracer maps of Sub- 
section 9.2. The region near the test section axis is considerably 
richer than the overall <J), as noted above. The structure of the 
rich swirler jets is not nearly as well-defined or symmetric. 

However, the hot and cold-flow data are roughly similar in a qualitative 
sense. 

In terms of evaluating the module wake enrichment factor F^/F q , 

a precise value cannot be obtained without weighting all the points 

in the wake according to their relative importance to the formation 

of NO^. This would involve detailed measurements not only of <J), 

but also the relative residence time at each point. The NO emissions 

x 

data matching of Subsection 3.2 produced a value of 2.5, For the 

contour maps of Figure 35, this corresponds to <p ^ 1.0. Generally 

it is seen that 0.7 £ <p £ 1.5. There are extensive regions of 

Figure 35b where <p is near unity. Upstream of these regions 

(Figure 35a) <J> is higher, and downstream (Figure 35c) (p is lower. 

Thus, in a rough sense, support is given to the approximation 

F /F £ 2.5. 

P ° 

Contour maps of mixing parameter (for the same axial stations as 
the <p maps) are shown in Figure 36, For all of the s maps the 
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contour Interval is 0.1 s unit. The shaded region in the upper left 
of each of these regions is the contour band 0.2<s<0.3. In 
Figure 36a there are patches where s is slightly less than 0.2, a 
couple of peaks which are suspected to be bad data points, then a 
sharp increase of s with radius to calculated values well above 
unity, which appear as random shading. These high values are not 
considered meaningful since the assumptions of Subsection 10.5 will 
not be valid in these regions. In regions where pure air is found, 
s will be undefined because the local fuel fraction is zero; for the 
purposes of Figure 36, s has been constrained to be zero at the test 
section wall. Figure 36b shows a more constant value of s at the 
upper left, with a more gradual increase to s»l. In Figure 36c, 
the s values are quite well behaved, increasing with radius toward 
three or four peaks where 0.8 < s < 1.0; beyond this, s decreases 
toward the s=0 constraint at the wall. These peaks seem to be 
vestiges of the swirler jet structure, information which was not 
visible in the corresponding <{> map (Figure 35c) . 

As in the case of the parameter F^/F^ it is not possible to 
obtain a precise value for the appropriate average module-wake mixing 
parameter s^. The NO^ emissions data matching of Subsection 3.2 
produced a value of 0.4. Generally it is seen that 0.25 £ s < 1+. 
There is a wide region where 0.3 < s < 0.5 in Figure 36c, although it 
is weighted toward a smaller test section radius than the regions of 
Figure 35 where $ ^ 1. Upstream, in Figures 36 b and a, the region 
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s ^ 0.4 is narrower. Thus, the s maps tend to support the approximation 

s 0.4. 

P * 

11. Summary and Conclusions 

This report has presented a model for predicting N0 X emissions from 
gas turbine combustors. The structure of the model is modular, and is 
not restricted to particular combustor geometries. It consists of a stirred 
reactor followed by plug flow reactors to which dilution air can be added. 

A modified Zeldovich scheme is used for the NO formation kinetics. Fuel- 
air nonuniformities are incorporated. Application of this model to a tubular 
premixed combustor, and the NASA swirl can combustor produced quite satis- 
factory agreement with available emissions measurements. Numerical calcula- 
tions of streamlines and fuel evaporation backed up the application of the 
model to the NASA swirl-can combustor. Visualization experiments were also 
conducted for better understanding. Confirming quantitative data came from 
experiments with cold-flow tracer and propane-burning gas sampling. A general 
model for NO^ formation such as the one described herein, supported by 
theoretical calculations and experiments, can be useful for the interpolation 
and extrapolation of test data, for examining the effects of design trade-offs, 
and for increasing the understanding of processes inside a gas turbine combustor. 

The major general conclusions are: 

(1) A modified Zeldovich kinetic scheme using equilibrium 
combustion products is sufficient for the calculation 
of N0 X formation under typical gas turbine conditions. 
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(2) The type of general N0 x model presented herein can 
extract information about entrainment and mixing in 
the combustor from engine emission data. 

(3) Low N0 x emissions require the combustion zone to be 
lean and very well mixed. 

(4) A measure of the mixedness in a combustor can be derived 
from time-averaged internal gas sample measurements. 

The following are the major conclusions specific to the NASA 
swirl-can combustor: 

(5) At the operating conditions of interest, combustion 

in the module wakes is only slightly lean and not very 
well mixed, but the residence time is reasonably short. 

(6) Fuel evaporation times are much shorter than the gas 
residence times and "droplet combustion" does not occur, 

(7) Low N0 x emission can be achieved only by making the 
module wakes leaner and considerably better mixed. 


I 
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